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Abstract 

Two-dimensional system of the fermions with the indirect Einstein phonon-exchange 
attraction and added local four-fermion interaction is considered. It is shown that in 
such a system at resulting attraction between particles a new nonsuperconducting phase 
arises along with the normal and superconducting phases. In this, called "abnormal 
normal", or pseudogap, phase the absolute value of the order parameter is finite but its 
phase is a random quantity. It is important that the new phase really exists at low car- 
rier density only, i.e. it shrinks with doping increasing in the case of phonon attraction. 
The relevance of the results for high-temperature superconductors is speculated. 
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1 Introduction 



The problem of an adequate description of the physical properties of high-temperature su- 
perconductors (HTSCs) still remains one of the actual problems of the modern solid state 
physics. It is connected with some peculiar properties of HTSCs. Among them there are 
such as quasi-2D character of electronic (and magnetic) properties, a relatively low and 
changeable carrier density rif and its influence on properties of HTSCs (see, for example, 
review 

Nowadays, one of the widely discussed questions on HTSCs is the problem of the so called 
pseudogap (or spin gap if magnetic subsystem of HTSCs is taken into account) 0, ^, |j, which 
is usually experimentally observed as a loss in the spectral weight of quasiparticle (or spin) 
excitations in normal state samples with lowered carrier density rif 0, §|, 0. Corresponding 
underdoped samples reveal some specific spectral, magnetic and thermodynamic pecularities 
which still continue to be not sufficiantly understood now. Moreover, the striking difference 
between the low (underdoped) and high (overdoped) density regions in HTSCs is increasingly 
debated and is considered as one of the very central and key questions in physics of the 
cuprates || || . 

The possibility of experimental changing of rif value in HTSCs puts a rather general 
theoretical problem of the description of the crossover from composite boson superfluidity 
(low rif) to Cooper pairing (large nf) when nf increases (in other words, a description of the 
transition from the so called underdoped regime to the overdoped one). Such a crossover 
was already studied for 3D and quasi-2D systems (see reviews [pi], [llj)- 2D case has been 



considered for the present at temperature T = only |T(], [12| what is connected with the 
Hohenberg-Mermin- Wagner theorem which forbids any homogeneous (long-range) order in 
pure 2D systems at T ^ due to the long-wave fluctuations of the charged order parameter 
(OP). 

The problem of the inhomogeneous condensate (Berezinskii-Kosterlitz-Thouless, or BKT, 
phase) formation was also considered despite of some difficulties in 2+1 relativistic field 
models [1131 where the fermion concentration effects are irrelevant. At the same time these 
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effects were studied in nonrelativistic model in |TJ|], for example, without taking into account 
the existence of the neutral OP p. Its consideration proves to be very important (see [15]]) and 
results in the formation of a separate equilibrium phase with p ^ which is located on the 
phase diagram of a system between normal and superconducting (here - BKT) ones. Due to 
fluctuations of the OP phase this new state of a system is of course also non-superconducting. 

In this paper an attempt is made to study the crossover as well as the above mentioned 
new phase formation possibility in 2D fermion system with both a more realistic indirect 
(phonon) and also a direct (local) four-fermion (4F-) interactions. Thus, the work is to a cer- 
tain extent a specific and non-trivial generalization of the preliminary short communication 



15| where this non-superconducting phase appearance was studied for 4F-case only and of 



the paper |16[ where Frohlich model was used for the investigation of the crossover at T = 0. 
As it will be seen in the boson exchange model (in contrast to the pure 4F-case), the new 
phase really exists when nf is rather small what allows to relate this result to underdoped 
HTSC compounds. But actually it is interesting to take into account a more real situation 
with an indirect attraction and some kind of local repulsion which may correspond to the 
short-range (screened) Coulomb interaction between carriers. In general case we, however, 
suppose that 4F-interaction can be repulsive and attractive as well. Besides, the case of total 
repulsion allows to explore the fermion- ant if ermion (electron-hole) pairing channel which in 
spite physical difference can be formally described by the same manner. 



2 Model and main equations 

Let us choose the simplest Hamiltonian density in the form: 

H{x) = -ipt(x) ^ - fij Mx) + H ph ((p(x)) + 
9 P h^l{x)i) a (x)^(x) - g 4 Fip\(x)^\(x)iJ l (x)^(x), (x = r,t), (1) 

where ip a {x) is a fermionic field with an effective mass m and spin a =f, J.; (i is the chemical 
potential of the fermions which fixes nf\ <p{x) is a phonon field operator, g p h and g^p are 
the electron-phonon and the 4F-interaction coupling constants, respectively. As we said 
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above, g^p can be positive (fermion-fermion attraction) or negative (fermion-antifermion 
attraction); in (jlj) we set Ti — kg — 1. 

In (|I|) H p h is the Hamiltonian of free phonons which can be described by the propagator 



where Q n = 2nnT (n is an integer) is the Matsubara frequency ||17|| . As it follows from 
(§), the propagator D(iQ n ) has been chosen in the simplest form with u being the Ein- 
stein (dispersionless) phonon frequency. It was done because of several reasons: first, this 
propagator gives a possibility to integrate the equations obtained; second, it is precisely 
the optic phonon and quadrupolar exciton modes with their relatively weak dispersion are 
widely considered as exchange bosons that can contribute into the hole-hole attraction in 
HTSCs [[[], [L8|, [19| , and third, the qualitative results concerning retardation effects do not 
strongly depend on the model studied. But on the other hand, the propagator (|2|) for the 
model under consideration can hardly be used for quantitative description of the cuprates 
and also their spin-wave branches which, as it is well-known, have linear dispersion law. 

It is important that the Hamiltonian ([!]) is invariant under symmetry transformations of 
two types |20|| , namely: 



z/v(z) - Mxy a(x) , 1>Ux) - i>t(x)e- ia W (3) 

and 

4)\{x) -> ip\{x)e- ia{x \ ip\{x) -> 4)\{x)e ia{ - x) (4) 

which must be taken into account. The phase a(x) in (||) and (|4|) is real. 

With the purpose to calculate the phase diagram of the system it is necessary to find 
its thermodynamic potential. It can be calculated by making use of the auxiliary bilocal 



field method (see, for example, |21|), which is a generalization of the standard Hubbard- 
Stratonovich one for the boson-exchange case. Then the grand partition function Z can be 
expressed through a path integral over the fermionic ip a {x) and the complex auxiliary fields 



(for example, 0(x, x') ~ < ip^(x)ip\(x') >). 



In the case of the model (|1|) it is convenient following the Ref . [ 52J to introduce the bispinor 

&(x) = (ip\(x), il>\(x), fy(x), ^(x)) (5) 



and its hermitian conjugate one which here are the analogous of the Nambu spinors |[23|| . 
After the substitution of ([Sp in (JTj) the Hamiltonian takes the form: 

A 2 

H(x) = -W(x)(— + f 1 ) 1 ® T z^>{x) - g ph W(x)I ® T z ^(x)(p(x) - 

g 4F ^(x)I <g> r z ^(x)^(x)I <g> r 2 ^(a;) + <^(a;)D _1 (x)^(a;), (6) 

where I ®t z is the direct product of the unit / and Pauli t z 2 <g> 2 matrices; is defined 
by (H). In such a representation of the Hamiltonian (ffl) and the field variables (El) the Fein- 



man diagram technique becomes applicable in its ordinary form ||22|| . Thus, after standard 
excluding of the boson field tp(x), the Lagrangian of the system can be expressed by the 
formula: 

A 2 

L(x u yt, x 2 , y 2 ) = W(x)[-d r + (— + n)I ® r z ]^(x) - 

® r 2 AT(xi, yi ; x 2 , y 2 )* (a*)*^)/ <g> r z . (7) 

The kernel .fT is the effective non-local inter-particle interaction function and will be explicitly 
defined in the momentum space below. 

In order to explore the pairing possibility in the system let us introduce the bilocal 
auxiliary field, or OP, 

<f>(xi,yi) = T z K(x 1 ,y 1 ] x 2 , y 2 )fy(x 2 )fy\y 2 ) = il ® T y $ ch (x u y x ) + t x ® r^^^i, y x ) (8) 

(the integration over x 2 and y 2 is assumed). Here 4>ch ~< wVm > an d 4>ms ~< V'jV'T > 
are electron-electron (charged) and electron-hole (insulating) OP, respectively (we neglect 
non-zero spin pairing). The auxiliary fields <p c h and <fii ns are responsible for the dynamical 
symmetry breaking (in according with (^) and (|4]), correspondingly). 
Adding to (|7p a zero term 

-[(j)(xi, y x ) - K(x x ,y X ] ^i,y[)^{x[) i if i {y' 1 )]K- 1 {x 1 ,y 1 ] x 2 , y 2 )[4>(x 2 , y 2 )- 



K(x 2 ,y 2 ;x , 2 ,y' 2 Mx' 2 )^(y , 2 )} 
to cancel the 4F-interaction, one could obtain the Lagrangian in the form: 



L{x 1 ,y 1 ;x 2 ,y 2 ) = ^{x^-dr + (- Vfj)I®r z 

2m 



{x 1 ,y 1 )]*{y 1 ) + 



(x 1 ,y 1 )K 1 (x 1 ,y 1 ;x 2 ,y 2 )(f)(x 2 ,y 2 ), 



Let us transform the expression for the kernel K; then in the momentum space it is 



(9) 



d 5 P(fp 1 d^p 2 

— K P ( qi ;q 2 )exp 



K(x 1 ,y 1 ;x 2 ,y 2 ) 



,xx + yt x 2 + y 2 . 
— — ) - ipi{xi -yi) - ip 2 {x 2 - y 2 ) 



2 2 

where pi = (pi,u>i)(i = 1,2) and P = [P,uj) designate the relative and the centre of mass 
momenta, respectively. By the definition the kernel Kp(pi,p 2 ) is in fact independent of P 
(so we omit index P henceforth) and aquires the simple form 



K(pi,p 2 ) = gl h D{p 1 - p 2 ) - g iF 



(10) 



which will be used in (f|) . The last expression evidently demonstrates that the total character 
of the effective inter-particle interaction as it always takes place in such a situation |24], p3| 
defined by the possible competition between the first (retarded) and second (retardless) 
terms in (|10D, or their common action. 

The partition function can be written as: 

Z = J ] VmV(j)V(j)* exp -f3 J L(^\ty,(f)*,(f))dxdy 



where Q[Q] is the thermodynamic potential which in the "leading order" is 

/3Q[g\ = -Tr Lng- 1 + ^Tr(0irV) 



where Tr includes 2D spatial r and "time" < r < (3 integrations as well as the standard 
trace operation. The full Green function of a system is 



g- 1 = -d T + r z + fi J I ® r z + 



(12) 



From (|TT| ) and flT2|) we arrive to the ^-equation (the Schwinger- Dyson one): 

——K(p;k,co)g(k,cu) = 0. (13) 
Substituting ([y|) into ([TTD one can obtain the expression for 



/3Q(S) = -TrLngT 1 + ^TiQKQ, 

The last is the well-known Cornwell-Jackiw-Tomboulis formula for the effective action in the 
one-loop approximation [25]. Using (|T3|) we can rewrite this expression in the form 

pn{Q) = -Tr[Ln£ + i^" 1 - 1]]. (14) 

As it was shown by Thouless et al. P6) (see also [|HJ) in 2D case it is natural to pass to 
a new parametrization of the OP (§) - its absolute value and the phase, namely: 

<Pch(x } y) = p ch (x, y)exp[-i(9(x) + %))], 

4>ins{x, y) = p ins (x, y)exp[-i(9(x) + 9(y))}, (15) 

where p c h and p ins are real. 

As it will be shown below, with the given kernel (Tjl]) there can arise only one (<j) c h or </>j ns ) 
OP. Therefore, it is necessary to make, simultaneously with (|T5j), the spinor transformation 
(in according with (|3]) and (|])) 

&(x) = x ] {x)exp{i6(x)I ® t z ), (16) 

&(x) = x ] {x)exp{i9{x)I <g> t s ), (17) 

(the spinor x( x ) is rea l an d formally corresponds to chargeless fermions). It is easy to see from 
(0), ( |T6"D and ( JT7|) that the phase dependences of the charged and insulating OPs are similar. 
Below we shall obtain ^-corrections for the (j) c h case, because the final equations for (pi ns will 
be the same up to substitution p c h — > Pi ns - The reason is that when K(pi,p 2 ) describes 
the attraction (charge pairing channel) the symmetry of the Lagrangian under operations 
(H) proves to be crucial for the representation (|16|); while when K(p 1 ,p 2 ) corresponds to 
the repulsion (chargeless, or insulating, pairing channel) the symmetry ([|) becomes already 
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important and the representation ( |TTD must be used as a "working" one. With this difference 
the rest of the calculations are almost identical and so we shall consider in detail the charge 
channel which is most interesting for metallic (superconducting) systems. 
In variables (|T6|) the Green function (|T2| ) transformes to 

/V 2 ^ 



g- 1 



-d T + I 



I®T z \d T + 



y2m 
Vfl 2 
2m 



+ \X + UTyP ch + 



'w 2 e v#v\ 



2m 



+ 



m 



J = G-\p ch ) 



E{d6). (18) 



Then using ( |l8l) supposing 6 gradients are small (the hydrodynamic approximation) and 
taking them into account up to the second order the effective potential (|14]) can be naturally 
divided it two parts: Q = Qkin{pch, V0) + Q po t{Pch) where in (V0) 2 approximation 



P^kiniPch, V0) = Tr 



1 



iGoSG S + r x ® J^p^GE + G£G£) 



(19) 



Assuming now that p c h(x, y) is homogeneous [] after somewhat tedious but otherwise straight- 
forward calculation one can obtain from (fO^): 



JWPcfc, A0) = g / dT / d M^^Pch(/i,T))(V0) i 



(20) 



where 



J(p,T, p c/l (/i,T)) 



1 

2tt' 



T 

7T 



Pch 



5 



4T2 0(^/47*) 



exp 



P 2 + /4 



T 

a; + p/2T 



^/ 2T cosh 



2 'x 2 



P 2 ,/4T 2 



(21) 



plays the role of neutral OP stiffness. Note that in comparison with the retardation free 4F- 



model [15j the last expression contains one more term, namely: the term with the derivative. 

The equation for the temperature Tbkt of the BKT transition can be written down after 
direct comparison of the kinetic term ( p0|) in the effective action with the Hamiltonian of 



1 Equations for p c h and pi ns will be obtained below and, as was shown in [ p^[ , it is an admissible approx- 
imation to put in them the value p c h (and pi ns ) independent of spatial and time variables. 



the 2D XY-model which has the formally identical form ||27|| . Hence it easy to conclude that 



/ Potion) ] T g f d 2 k 



m=— oo 



(2tt) 



TC 

— J(p,T BKT , p ch (p,T BKT )) = T BKT . (22) 

The essential difference of this equation obtained from that for the XY-model is its sponta- 
neous dependence on p (or rif) and p c h- 

To complete the set of self-consitent equations which allow to trace an explicit dependence 
of T B kt on rif, the equations for p c h and p have also to be given. In particular, the equation 
for p c h{iuj Tl ) is nothing else but ( |13"1) with V6* = 0, i.e. the Green function G of the neutral 
fermions substitutes Q, so that ( |13"1) in frequency-momentum represantation takes the form 

-pchiiuJm) \ K{u} n ,U} m ) 

+Pins{iuJ m ) ) U ™ + + Pch( iuJ "i) + Pinsi^m) ' 

(23) 

where uj n = (2n + l)nT is the Matsubara fermionic frequency and kernel K{uj m ,u n ) is 
defined above. We cited the final equations for both OPs, p c h and p ins , in order to show only 
that they indeed are the same but alternative if the kernel K changes the sign. 

An analitical solution of these equations, as well as obtaining both the equation (|22]) and 
the number equation needed is only possible if one supposes that p c h{iuj n ) does not depend 
on the Matsubara frequencies (see footnote on p. 8.). 

Making use of this approximation the number equation which follows from the ordinary 
condition V~ 1 dQ[Q]/dp = —rif (V is a volume of a system) and is crucial for crossover 
description has to be added to (1221) and (E3p for self-consistency; so one comes to 



p 2 + p 2 ch + p + 2T\n 



exp 



P 2 + Plh 



T 



2e F , (24) 



where e F = irnf/m is the Fermi energy of free 2D fermions with the simplest quadratic 
dispersion law. Thus, in the case under consideration all unknown quantities p c h, p and 
Tbkt are the explicit functions of rif. 
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3 Analysis of the solutions 

Unlike the usual (with T-independed unit vector) XY-model, in the superconducting one 
there exist two critical temperatures: T p where formally the complete OP given by (|8]) arises 
but its phase is a random quantity, i.e. < (f)(x,y) >= 0, and another one, Tbkt < T p , 
where the phase of the OP becomes ordered, so that < <j)(x,y) >^ 0. In other words, 
the temperature T p is in fact the temperature of appearance of neutral OP only which has 
discrete symmetry and thus is not at variance with general theorems. Recall, that according 
to the equations obtained above both these temperatures, what is important, directly depend 
on the carrier density in the system. 

Critical temperature T p can be found from (|2"T|)-(|2~3|) by putting p ch = (what in ac- 
cordance with derivation of these equations corresponds, in fact, to the mean- field approx- 
imation). As a result, a 2D metal with temperature decreasing passes from normal phase 
(T > T p ) to another one where averaged homogeneous (charged) OP < (f)(x,y) >= 0, or, 
what is the same, superconductivity is absent, but chargeless OP p c h ^ 0. It is evident 
that the pseudogap is formed just in the temperature region Tbkt < T < T p , because, as 
follows from the above formulas (see, for instance, (^l|)- (|2^) ), p ch = p c h(T) enters all spec- 
tral characteristics of 2D metal in the same way as the superconducting gap A(T) enters 
into corresponding expressions for ordinary superconductors. It justifies why this new phase 
can be called the "abnormal normal" phase or better pseudogap one. The density of states 
near ep in the pseudogap is definitely less than in the normal phase, but is not to be equal 
zero as in superconducting one. The latter has to be checked by direct calculation of the 
one-particle fermion Green function what is most likely a separate problem which we do not 
touched upon here. 

The phase diagram of the system can be found from the equations (|2l| ) - (|2~4]) . There 
are different behaviours of T p (rif) and TsKT^f) for various correlation between interaction 
constants. 

2 Because of p c h and pi na can not exist simultaneously (see (p3|)) the index p means the only OP, which 
appears at finite sign of total interaction. 
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1) g^F > 0,g p h = (retardless interaction). 

This case has been partly analyzed in Ref . . It corresponds to fermion-fermion pairing 
due to the local interaction. Note (see (|23|)) that in the case of resulting attraction be- 
tween fermions fermion-antifermion (or electron- hole) insulating pairing channel is absent, 
i.e. pi ns = 0. The phase diagram for this case is presented in fig. 1. It shows that the 
abnormal normal phase exists at any concentration value n/ and the temperature width of 
this phase region weakly increases with rif increasing, and BKT phase always begin to form 
when p ch (T BKT ) is finite. 

At ep — > the temperature of BKT phase formation is defined by equality T BKT = £f/2, 
and T p as function of rif can be found from the equation T p \n(T p / 'ep) = W exp(— 2/g^pm), 
which follows from (|23| ) (W is the conduction band width). 

2) giF > 0,g ph ^ 0. 

The situation here is almost the same as previous one. The presence of the indirect 
interaction leads to the effective growth of the effective 4F-interaction constant g\ F > g^p 
(see (^)). The latter in one's turn simply results in increasing of the region between T p and 
Tbkt (fig-2) and keeps the form shown on fig.l. 

3) <?4F = 0,g p h 7^ (a pure indirect interaction). 

This is one of the most interesting cases because it corresponds to the widely accepted 
electron-phonon (or BCS-Bogolyubov-Eliashberg) model of superconductivity. The numeri- 
cal calculations of the phase diagram is presented in fig. 2. This diagram shows that compar- 
atively large region with the abnormal normal (pseudogap) phase exists at rather low carrier 
concentrations only and its temperature area shrinks at rif —>■ oo. Such a behaviour qualita- 
tively agrees with that one which takes place in real HTSCs samples || ||, [7|, || demonstrating 
that pseudogap (and spin gap also) region is observed in underdoped samples. 

Indeed, it is not difficult to make certain that the assymptotics for T p (n/) and Tb^t(^/) 
have the following forms: 

i) when ratio ep/u 1 (very low fermion density, or local pair case) the first of 
them satisfies the equation T p \n(T p / ep) = u exp(— An / g^ h m) which immediatedly results 
in dT p (nf)/drif\ nf ^o — > oo. At the same time the temperature Tbkt at n>f — ► has another 
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carrier density dependence and as above Tbkt = cf/2 what simply means that here it again 
is equal to the number of composite bosons; in this density region T p /Tbkt S> 1 (such an 
inequality is also correct for the pure 4F-interaction) . 

ii) in the case ep/ojQ ^> 1 (very large fermion density, or Cooper pair case) one easily 
arrives to the standard BCS value: T p = {2^ujq/t[) exp(— 2ix / g 2 h m) = T^q S = (2j/tt)Abcs 
(Abcs is the usual one-particle BCS gap at T = 0). In other words, the temperature T p in 
this limit becomes equal to its BCS value 0. The T BKT asymptotics for this case is not so 
evident and requires more detailed consideration. 

First of all, it is naturally to suppose that for large rif value Tbkt — ► T p . Then it is 
necessary to check the dependence of p on T as T — ► T p . For that the equation (|23| ) can be 
transformed to: 



9%™ Jo V J x 2 + p 2 h/AT 2 2(Jx 2 + p 2 ch /AT 2 -uj /2T) 



2tt r°° I tanh ^/a; 2 + p 2 h /4T 2 tanh ^x 2 + p 2 ch /4T 2 - tanh(w /2T) 
dx 



tanh ^x 2 + p 2 ch /AT 2 + tanh(o; /2T) \ 
2^x 2 + p 2 ch /AT 2 + u Q /2T) ) 
(where it was used that in this concentration region, the ratio p/2T p ~ ep/2T p ^> 1 because 
of//~e F 0,[TJ,|g,[r6|). 

On account of usually u>o/2T p ^> 1, only very small x give the main contribution to the 
integral ([25|) (it is seen from the limit p/2T p —>■ when e^/cjo — » oo). Therefore the latter 
expression takes the approximate form: 



g 2 h m Jo y ^ x 2 + plJAT 2 x + u /2T ) ' 
On the other hand, the condition p c h = in (p6|) leads to the equation 



2n r°° , / tanhx 1 \ , . 

dx — — - lnm ■ (27) 



g; h m Jo \ x x + uj /2T Pchi 



3 Being equal (in mean field approximation only) these temperatures (T p and T^Jg) are in fact different: if 
Trcs immediately falls down to zero as fluctuations are taken into account, T p does not and is renormalized 
only. 
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for T p . From (|26| ) and ( p?| ) it directly follows that 



°° / tanh x tanh ^x 2 + pg fe /4T 2 

tXX' I 
\ .T 



* 2 + f?jATZ 



T 
ln^. 

T 



Then using the approximation 



tanh Jx 2 + /4/4T 2 



1-3-V + /4/4T 2 ], x<l; 



/a: 2 + p 2 ,/4T 2 | x" 1 - p c \/8TV, 
one directly comes to the expression needed: 



x > 1, 



p ch (T) ~ 2.Q2T p JT p /T - 1. 



(28) 



Recall that the generally accepted 3D result is A B cs(T) = 3.0QT B( f s y T B( f s /T — 1 |HJ and 
this small difference can be explained by the above approximation what, however, is suitable 
for the following below qualitative discussion (see next Section). 



The dependence ( |2"8"| ) has to be substituted in equation (|2"2"|). And again because of 
fx/2T B KT — ^f/^Tbkt > 1 and p c h(T BKT )/2T BK T "C 1 when T B kt — ► T p this equation can 
be written as 



e F 



AT, 



BKT 



Plh 







d(p 2 ch /4T 2 BKT ) 



dx 



( 1 



V cosh 2 



x 



cosh 2 + p 2 ,/4T 2 KT 
At last, using expansion in p c ^ / 2T B xt in integral fl2"9|), the latter can be transformed to 



(29) 



a^F ( Pch V 
^T B kt \^T B ktJ 



(30) 



where the numerical constant 



a = I dx 
h 



tanh 2 x — x 1 tanh x + 1 



1.98. 



2x 2 cosh x 

Combining now ( ^8|) and (|30"D one comes to the final simple relation between T Pch and T B kt 
for the large carrier density: 



T BKT ~TJ1- 1.17 JT p /e F 
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i.e. Tbkt as a function of n/ really approaches T p (or T^f s ) (see fig. 2). 

As to crossover region defined by the equality fi ~ it is easy to convince from the 
equations (|21"|)-(|24|) and fig. 2 that it corresponds to the densities when the temperatures 
T p are essentially different. It is important that because of relatively low for the phonon 
case value of the bound pair-state energy and so very small region of negative rif Hl6| , the 
behaviour Tsxrinf) ~ €f hardly corresponds to Bose-Einstein condensation and in fact 
takes place at \i > 0. 

4) 94F < 0,9 P h ^ 0, but g 2 ph » \g 4F \. 

This condition provides the total fermion-fermion attraction channel only, so p c h ^ 
(see below) . As it was said above the local 4F-repulsion qualitatively can correspond to the 
screened Coulomb repulsion. In this situation the cut parameter must be introduced to avoid 
the divergence in (p3|). It is interesting to consider two situations: i) the cut parameter (the 
boundary Coulomb frecuency u c ) goes to infinity (i.e. small concentration or local pairing) 
and ii) when uj c is large but finite u c » u but g ph » |<74f|. Probably it is the situation 
what is intimately related to the real HTSCs (and superconductors at all). 

i) The case 1) is restored in general but the effective 4F-constant g\ F < 0. It is important 
that Tbkt preserves its linear assymptotics at small rif. 

ii) In this case g^p in (^) can be substituted by g±FD{iu) n ) with ui = oj c (the propagator 
D is defined by ([?])). The situation is similar to the case 3) but the effective coupling constant 
is smaller. Such a decreasing leads to the narrowing of the abnormal normal phase region 
because of lowering of the temperature T p . In particular, it is not difficult to obtain the well 
known Tolmachev logarithmic correction to T p (see, for example, f22|1 ): 

T p = ^^-exp{ — 2 1 



where fi* = g^N^O)/ (1 + g^N^O) In^ep/ujo)) (iV(0) is the density of states at the Fermi 
surface) . 
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4 Conclusion 



The model proposed to describe the possible two-stage superconducting phase transition in 
2D (and quasi-2D) metallic systems is in fact very simplified in order to investigate their 
most typical and general features. All the more surprisingly that it catches some essential 
details which are characteristic for underdoped HTSC copper oxides. In particular, the 
experimental data demonstrates |28], ^] that i) indeed for low rif the critical temperature 
T c is proportional to rif (what is simply ep), ii) T c shows saturation when rif approaches 
so called "optimal doping" (i.e. carrier concentration when T c as function of rif reaches its 
highest possible in given compound value), iii) the ratio T c /ep in these and other "exotic" 
superconductors is as high as 10~ 2 — 10 _1 what independently points out on rather small 
Fermi energy, etc (for details see p^ ). 

One would think that the pecularities mentioned receive their natural explanation on 
the basis of the model of metal with indirect inter-fermion interaction if the temperature 
Tbkt is implied as critical one T c (this is justified for pure 2D systems [[30]). In quasi-2D 
model because of the third spatial direction and the phase fluctuation stabilization the true 
temperature of ordinary homogeneous ordering arises |29|, |3l| (see also 

As regards the second temperature, T p , it usually introduced by empirically as some 
temperature point T* where observable spectral (or magnetic) properties of HTSC begin to 
deviate appreciably from their standard for normal metallic state behaviour || [], [7], || . As 
a rule such a deviation is connected with appearance of fluctuating (short- living) pairs. We, 
however, showed that some finite number of these pairs does exist and is to formed at definite 
temperature due to phase transition between normal and pseudogap (also normal) phases. 
The only difference with supposed dependence T* on the density of doped holes consists of 
assymptotics at rif decreasing: we have obtained that this temperature is also reduced while 
usually (see, for example, [p9jD T* is presumably depicted as such one that incriases with 
rif decreases. It seems that the latter has no satisfactory grounds. Nevertheless it must 
be stressed that the above limit T p (rif) — > when rif — > can not be also considered as 
sufficiently regular because of the growth here of the neutral OP fluctuations contributions 
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of which was not taken into account and which become rather important at small rif. 

At last the model under consideration qualitatively correctly describes the explicit nar- 
rowing of the pseudogap phase area at carrier density increase (such a diminution results 
in rather rapid reapproachment the temperature T c (Tbkt) and T p and their experimental 
confluence (indistinguishableness) in BCS limit. 

Some important problems remain open and are to be solved. Among them there are: 
more complete and deep development of the model which has to consider different kinds 
of the dispersion laws for the intermediate bosons; more careful taking into account of the 
Coulombic repulsion; neutral OP fluctuations, especially at low n^; generalization of the 
approach on the case of non-isotopic pairing. On the other hand, high-T c compounds must 
be investigated in the frame of more realistic model that such their pecularities as magnetism 
of cuprate layers, non-qudratic free carrier dispersion law with van Hove singularities in 
the hole density of states. One of the most interesting problem is to obtain doping and 
temperature effective action which is equivalent to Ginzburg-Landau potential because in 
many cases the phenomenology is more preferable. 

We would like to thank Prof. V.P. Gusynin, Drs. S.G. Sharapov and LA. Shovkovy for 
extremely interesting and useful debates concerning the questions rised above. Expecially 
we are indebted LA. Shovkovy for reading the manuscript and comments. 
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Figures caption 



Fig.l 

Phase T - rif diagram of 2D metal with 4F fermion attraction. The lines correspond to the 
functions T p (rif) (the upper line) and T B KT(n f ) (the lower one) at A = 0.5. The figures I, II 
and III show the regions of the normal, abnormal normal (pseudogap) and superconducting 
phases, respectively. 
Fig.2 

Phase T - nf diagram of 2D metal with indirect inter-carrier attraction. Similarly to the 
Fig.l, the lines correspond to the functions T p {nf) and Tbkt{u s ) and separate the same 
regions (A = 0.5). 



19 



